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Abstract 


It is shown that in the framework of analytical confinement, when 
quark and gluon propagators are induced by an vacuum selfdual gluon 
field with constant strength, the masses of meson with quantum num¬ 
bers Q = and quark constituents mi, m 2 are described with rea¬ 
sonable accuracy by the formula 


Mq {mi , m2) = (mi -r m2) 


1 + 


A, 


Q 


{mf + 1.13mim2 -|- m^) 


2 ', 0.625 


where a constant positive parameter Aq is unique for all mesons with 
quantum numbers Q = . 

Sets of mesons = 0~, 1“, 0+, I"'', 2+, 3“ and different flavors 
constituent quarks {u = d, s, ,c,b) are considered. 


1 Introduction. 

Our basic point of view is that the conhnement and hadronization of qnarks 
take place on the same distances. It means that the knowledge of behavior of 
qnark and glnon propagators on these distances give us direct way to calculate 
meson masses using the Bethe-Salpeter equation in ladder approximation if 
the coupling constant is small enough. 
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Within the standard approaches of QCD the quark and gluon propagators 
in the conhnement region are not obtained yet, so that different symmetry 
arguments leading to some effective theory are used to avoid direct calcula¬ 
tions. 

In the papers [1, 3] we developed the approach based on the assumption 
that the QCD vacuum is realized by the self-dual homogeneous vacuum gluon 
held which is the classical solution of the Yang-Mills equations. The propa¬ 
gators of quarks and gluons in this held can be calculated in the explicit form 
and they are entire analytic functions in p^-complex plain. It means that the 
self-dual homogeneous vacuum gluon held leads to the analytic conhnement 
of quarks. Manifestation of this gluon conhguration in the spectrum and 
weak decays of light mesons, their exited states, glueballs were studied in 
papers [2, 4]. Our calculations showed that the self-dual homogeneous gluon 
held can be considered as a good candidate to realize QCD gluon vacuum. 

However any calculations in this model are quite cumbersome so that it 
is difficult to ” see forest among trees”. I always wanted to hnd some simple 
formulas which could have some general physical consequences. I think this 
idea is realized in present paper. 

Physical arguments are the following. Bound states are formed on large 
distances where conhnement plays the main role. So that the correct descrip¬ 
tion of conhnement should lead to reasonable description of bound states. 
Mesons as quark-antiquark bound states are characterized by diherent quan¬ 
tum numbers which should be taken into account for detailed description 
of meson characteristics. However if the analytical conhnement really takes 
place then one can get some general dynamic properties of meson spectrum 
which weakly depend on quantum numbers. 

We use the functional integral techniques to calculate mass spectrum of 
bound states. In this techniques representation of the Bethe-Salpeter kernel 
with propagators, which are entire analytical functions of the gaussian type, 
has some specihc features leading to consequences which can be checked on 
the real meson spectrum. Namely this kernel is product of a polynomial 
depending on quantum numbers and an exponent which is function of quark 
masses and mass of the bound state under consideration and which gives 
the main contribution into integral. Thus we can expect that in the lowest 
semiquantative approximation the meson masses are dehned by an univer¬ 
sal function depending on masses of constituent quarks and dependence on 
quantum numbers can be approximated by a constant. 
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Let us consider the Particle Data Group [6]. One can extract some sets 
of mesons having the same quantum numbers Q = J^, but different flavors 
constitutes. There are hve sets of mesons having as constituents all four 
quarks u = d, s, c, b. These sets are: pseudoscalar 0“ - ten mesons, vector 
1“ - eight mesons, scalar O’*" - four mesons, axial - five mesons, tensor 2~^ 
- six mesons (see Section 3). In addition we consider tensor multiplet 3“, 
containing three light mesons. 

Thus we expect that the masses of these mesons with quantum numbers 
Q = and quark constituents mi, m 2 are described by the formula 


MQ(mi,m2) = A^(mi,m2, Aq), (1) 

where Af (mi, m 2 , Aq) is an universal function, form of which does not depend 
on quantum state. A constant parameter Aq dehnes quantum state and it 
is unique for all mesons with quantum number Q = J^. 

The aim of this paper is to show that the representation (1) really takes 
place and to clarify additional conditions providing this formula. 


2 Bethe-Salpeter equation, 

analytical confinement and meson masses 

The equation on the mass M of a bound state of two constituent quarks with 
masses mi and m 2 in framework of the Bethe-Salpeter equation in ladder 
approximation looks as (see [1]) 

^ = dyidy2VQ{yi)IlQ{yi - y2\p)VQ{y2). ( 2 ) 

Here the polarization operator is 

AQiyi-y2\p) (3) 

= J dx e*(^^^Tr [rQS'i(a; -h p,2{yi - y 2 ))TQS 2 {x - pii^yi - 1/2))] 

where the quark-antiquark vertex Tq dehnes the quantum numbers of the 
bound state under consideration. 

All calculations will be performed in the Euclidean space R^. 
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The quark propagator is chosen in the form induced by a vacuum self-dual 
gluon held with constant strength (see details in [1]) 


S{y,m) 




^S{y,m,u) ■ e 


1-u 
1 + u 


(4) 


The parameter A dehnes the conhnement scale. The explicit form of the 
polynomial S{y,m,u) over variables y and m will be not important in our 
calculations here but it can be found in [1], 

The other parameters are 




M 


y = 


mi + m2' 


hi — 


mi 


mi + m2 


A 2 — 


m 2 


mi + m2 


The vertex function Vq is determined by the solution of the Bethe- 
Salpeter equation, but with acceptable accuracy it can be approximated on 
large distances ~ by 


VQ{y)r^D{y)r^V e-^y\ 


(5) 


The QCD coupling constant is dehned by the standard way 

.2 


a, = 


El 

dTT 


( 6 ) 


Thus the mass of a bound state with quantum number Q is dehned by 
the equation 

1 ^ J dy V^e~^ J dx e''‘'y'''^TT[TQSi{x + iJ, 2 y)TQS 2 {x - niy)] (7) 


Here Pq{ui,U 2 , y, yi, ^ 2 ) is a polynomial over parameters y, yi, y 2 and its 
explicit form is dehned by the spin structure of verteces and qnark propaga¬ 
tors. The qnantnm nnmbers Q = of bonnd states are dehned by this 
polynomials. The explicit form of this polynomials will be not important in 
our arguments. 
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The function 


( 8 ) 

_ 2 ^1^2 + 2(/iini + nlu2) f4 , + 

^ U 1 + U 2 + 2 2 VI —Ml/ 2 VI —^2/ 

plays the main role in the equation (7) and its behavior dehnes the meson 
mass spectrum. It is important that this function depends on masses of 
constituent quarks rui and m 2 and does not depend on quantum numbers of 
bound states, so that we can hope that the main features of a set of mesons 
with the same quantum number Q = are dehned by masses of constituent 
quarks only. 

Let us consider the behavior of the function E. Let us note bounds of 
changing of our parameters: 

M 

0</i=-<1 h- 2.5; /ii + /i2 = l. (9) 

mi + m2 

In the case u = —— < 1 (the meson mass is less then sum of masses 
of constituent quarks mi + m 2 ), the function E is negative and main contri¬ 
bution into the integral (7) gives the vicinity of the point Ui = 0, M 2 = 0. 
Then 

and the solution of the mass equation requires strong coupling regime > 1. 

In the case a = —^— > 1 (the meson mass M is more then sum of 
masses of constituent quarks mi + m 2 ), the function E is positive and has a 
positive maximum in a point 0 < m® < 1, 0 < w® < 1. Thus one can write 
approximately 


2 

JJ dUidU2PQ{Ui, U2, /i, /Ui, /i2)e 




where 


S{M,mi,m 2 ) = i7(/i,/ii,/i 2 ) = max i7(/i,/ii,/i 2 , Wi, M 2 )- (10) 
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The parameter Cq contains all information about quantum numbers of states 
under consideration. We can suppose that Cq depends weakly on mass pa¬ 
rameters /i, /ii, /i 2 - 

Calculations of the maximum in (10) in the region (9) give the approxi¬ 
mate formula 


(mi -F m2f£{M,mi,m2) 

K. [l.37mi -I- 1.55mim2 -l- 1.37m2 


M 


mi -|- m2 



Thus the equation (7) can be approximated by 


TT 


( 11 ) 


( 12 ) 


The solution of this mass equation requires weak coupling regime «* < 1. 
Finally the formula giving the mass of meson, which is in a quantum state 
Q = and is constructing by quarks with masses mi and m 2 , is 


Mq( mi, m 2 ) (mi m 2 ) 


1 + 


__ 

(mf -I- 1.13mim2 -l- ml)^'®^^ 


(13) 


Here the positive constant Aq is the same for all mesons with a given quantum 
number Q = J^. 

Now we should check this formula on real meson spectrum. 


3 Meson masses. 


Let us choose in the Table Particle Group Data [6] mesons with the same 
quantum numbers but with different flavors masses. Mesons can be rep¬ 
resented in the form 


Q = J^ 


/ uVqU 

uVqS 

uVqC 

uVqb 


sVqS 

sVqc 

sVob 



cVqc 

cVqI 

V 



hVqh 


(14) 


6 



We can choose only five sets of mesons with fixed having all four con¬ 
stituent quarks u = d, s, c, b 


0" 


1“ 


0 + 


1 + 


2 + 


/ r/(547) 77(494) D(1869) 5(5279) \ 

V(957) 5,(1968) 5,(5369) 

77e(2979) 5^(6400) ’ 

V Vb^OO) 

/ uj{782) K*{892) 5*(2007) 5*(5325) \ 

0(1020) 5! 2112 

J/0(31OO) - ’ 

V T(9460) 

f /o(980) - - - \ 

/o(1370) 

Xco(3415) 

V Xw( 9893) y 

/ ai(1260) i7i(1270H- 1400) - - \ 

/i(1420) 

Xci(3510) 

V Xw(9892) y 

/ 02(1270) 77* (1430) 5* (2460) - \ 

02(1525) 

Xc2(3556) 

V 762(9912) 


In addition we consider three mesons 3 : 

y ^3(1670) 773*(1780) \ 
03(1850) )■ 


We want to show that the masses of all these 36 mesons with reasonable 
accuracy are described the formula (13). Our parameters which should be 
determined by fitting are masses of constituent quarks m,, rric, rrib 

(4 parameters) and parameters Ap, Ay, As, A a, Ap, Ap (6 parameters). 
Our calculations consist of the following steps. 

The first step. We calculate the minimum 


On = min 


^Q,mf ,771 f {mf,m,,AQ) 


(Mq: 


.mf,mf yexp 


(15) 
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for each Q = P, V, S, A, D. These calculations give the parameters 
Ap, Ay, Asi Aa, Ap. On this step quark masses mj turned out to be 
a little bit different for different Q. 

The second step. The quark masses ruf are determined by calculation of 
the minimum 


= min 


E E 

Q q:u,s,c,b 


1 - 


MqQ{mq,mq,AQ) 


{M,q) 


exp 


( 16 ) 


with parameters Aq are hxed by (15). This calculation give the constituent 
quark masses nif = niu, nig, nic, mb- 

The third step. We calculate the parameter Ap 


Up = min 

Arr ^ 


1 - 


^T,mf ,m,p {mf,mf>,AQ) 

(Mt, 


rrif^m^i /exp 


(17) 


using the constituent quark masses ruf = rric, mb- 

The fourth step. We calculate the accuracy of our approximation for each 
given quantum number Q 


V mi,m2:w,5,c,o 


Mq ,mi ,m2 Aq) 


(Mq 


,m2 / exp 


( 18 ) 


where Nq is the number of known mesons with quantum number Q. 
The results are 


rrif 

Mev 

rriu 

260 

rus 

434 

rric 

1506 

rUb 

4732 



Q 

1 

o 

0 = 1- 

+ 

o 

/I = 1+ 

D = 2+ 

T = 3- 

AQ{Gev‘^) 

0.0216 

0.217 

0.249 

0.527 

0.618 

0.838 

Oq 

0.011 

0.0026 

0.012 

0.0023 

0.0019 

0.00016 


and are given in the Figures 1-6. 

Our results can be formulated in the form: 




• The formula (13) correctly describes the mass dependence of mesons 
on masses of constituent quarks. 

• The main features of meson spectrum are defined by the function 
T(M, mi,m 2 ) which does not depend on quantum numbers of a me¬ 
son state but on the quark constituent masses. 

• The coupling constant should be small enough < 1. 

• The mass of a meson should be more then sum of masses of constituent 
quarks M > mi + m 2 . 

• Masses of constituent quarks mq are close to the standard customary 
values. 

• This consideration is too rough for light pseudoscalar mesons tt, K, t], rj'. 

• There is large difference between parameters Aq. 

This paper is supported by RFBR grant N° 04 — 0217370. 
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Figure 1: Masses of pseudoscalar mesons P = 0~. Dotted line is M 
mi+ m 2 = m, solid line is M = M (jj-: T’ 
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Figure 2: Masses of vector mesons V = 1 . Dotted line is M = mi+m 2 = m, 
solid line is M = Af AyJ 
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Figure 3: Masses of scalar mesons S = O’*". Dotted line is M = mi +m 2 = m, 
solid line is M = Af (f, f, 2 I 5 ) 
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Figure 4: Masses of axial vector mesons A = 1+. Dotted line is M = 
mi + m 2 = m, solid line is M = Af Y’ 
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Figure 5: Masses of tensor mesons D = 2'^. Dotted line is M = mi+m 2 = m, 
solid line is M = M Y’ 
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Figure 6 : Masses of tensor mesons T = 3“. Dotted line is M = mi+m 2 = m, 
solid line is M = (f, f, 
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